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Abstract. We present the results of TNG near-IR low resolution spectroscopy of two (S7F5_254 and S7F5_45) 
sources belonging to a complete sample of 15 EROs with K’< 18 and R-K’> 5 selected from the MUNICS Survey. 
Both the spectra show a sharp drop in the continuum which can be ascribed only to the Balmer break. This 
places them at 1.2< t < 1.5. Their rest-frame 2 : = 1.2 K-band absolute magnitude is Mjc ~ —26.6 (L~ 7L*). The 
comparison of the spectra and the photometric data with a grid of synthetic template spectra provides a redshift 
2 ~ 1.22^0 2 ^ for S7F5_254 and 2 ~ 1.46 ± 0.02 for S7F5_45. The resulting lower limits to their stellar mass are 
•At stars d X 10“ Mq and MTtTrs = 4 X lO^*^ Mq. The minimum age of the last burst of star formation in 
S7F5_254 is 3.5 Gyr while it is 0.5 Gyr in S7F5_45 implying a minimum formation redshift Zf>3.5 and Zf>2 for 
the two EROs respectively. 
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1. Introduction 

The existence of massive (e.g. Mstars > lO^^M©) galax¬ 
ies at z > 1 and their spatial density provide crucial 
constraints on the current galaxy formation and evolu¬ 
tion models. The current ACDM hierarchical merging sce¬ 
nario predicts, indeed, that massive galaxies are assembled 
by means of subsequent mergers of disk galaxies largely 
occurring at z < 1.5 (e.g. Kauffmann 1996; Cole et al. 
2000; Baugh et al. 2002). In this scenario the higher the 
stellar mass of galaxies, the lower is the probability of 
their existence at increasing redshift. For instance, the 
predicted space density of Mstars > IO^^Mq galaxies is 
<10“®/i^Mpc“^ at z ^ 1, i.e. <0.02 galaxies arcmin“^ at 
1.2 < z < 1.5 {h = 0.7) (Baugh et al. 2002). An alter¬ 
native scenario is the “monolithic collapse” where, in the 
traditional form, even the most massive galaxies formed 
at z/>3 from the collapse of proto-galactic gas clouds. 
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Almost all of their stellar mass is produced in a single 
short episode of star formation. Their evolution should 
then follow a passive aging (Tinsley 1977; Bruzual and 
Kron 1980) as suggested by the observed properties of lo¬ 
cal spheroids (e.g. Renzini & Cimatti 1999; Peebles 2001; 
for recent reviews). A modern and more realistic view of 
the single-collapse model predicts that the star formation 
in elliptical galaxies lasts longer than 10® yr (e.g. Jimenez 
et al. 1999). In this scenario the spatial density of massive 
evolved galaxies is almost constant over a wide redshift 
range (0 < z < Zf). This is the main reason why in the 
last few years, the search for z > 1 massive galaxies has 
been attempted. Samples of candidates of evolved galax¬ 
ies at z > 1 have been selected on the basis of optical-IR 
colors (e.g. R-K> 5), the so called Extremely Red Objects 
(EROs; e.g. Thompson et al. 1999; Scodeggio & Silva 2000; 
Daddi et al. 2000, McCarthy et al. 2001; Martini 2001; 
Roche et al. 2002). Because of the age-dust degeneracy, 
both z > 1 old stellar systems and dusty star forming ob¬ 
jects spanning a wider range of redshift can populate the 
ERO regime in color space (Cimatti et al. 1999; Dey et 
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al. 1999; Mannucci et al. 2002). The recent estimates by 
Cimatti et al. (2002), based on optical spectroscopy of a 
complete sample of K< 20 galaxies over 50 arcmin^ and by 
Mannucci et al. (2002) based on the color-color diagram, 
show that the two classes are about equally populated and 
that a substantial population of passively evolved galax¬ 
ies at 2 ~ 1 exists (see also Franceschini et al. 1998). At 
higher redshift (z > 1), where the two paradigms of galaxy 
formation could be better distinguished, the existence of 
massive and/or evolved galaxies has not yet been proved. 
The detection of a massive disk galaxy at z = 1.3 by van 
Dokkum & Stanford (2001) could cast doubts on the cur¬ 
rent hierarchical models since they do not predict such 
regular, massive disk systems at z > 1. However, it could 
be an exceptional case. Evolved galaxies have been dis¬ 
covered at z > 1.5 (Soifer et al. 1999; Benitez et al. 1999; 
Stiavelli et al. 1999) but their small stellar masses do not 
provide severe constraints. 

The search for massive evolved galaxies at z > 1 needs 
both a wide angle K-band survey and spectroscopy in the 
near-IR. Indeed, the expected apparent magnitude of a 
■Mstars > IO^^Mq galaxy at 1 < z < 1.5 is K~ 18 — 18.5 
(e.g. Kauffmann & Chariot 1998). The observed surface 
density of K< 18 EROs is 0.05-0.1 arcmin”^ (e.g. Daddi 
et al. 2000) and about half of them can be dusty galax¬ 
ies (Cimatti et al. 2002). Thus, a reasonable probability 
to find z > 1 massive evolved galaxies can be achieved 
over areas larger than few hundreds arcmin^. Moreover, 
spectroscopy in the near-IR is needed to make feasible the 
detection of the 4000A break at these redshift. 

In this paper we report the analysis of two K’< 18.0 
un-lensed EROs spectroscopically observed in the near- 
IR as a part of the ongoing project TESIS (TNG EROs 
Spectroscopic Identification Survey; Saracco et al. 2002), 
aimed at searching for Mstars > IO^^Mq evolved galaxies 
at z > 1. Throughout this paper we assume Ho=65 km 
s“^ Mpc“^, Ho = 0-3 and Aq = 0.7. 

2. Observations and data reduction 

The two EROs presented here belong to a complete sam¬ 
ple of 15 candidates massive z > 1 evolved galaxies se¬ 
lected from two no-adjacent areas of about 180 arcmin^ 
each from the Munich Near-IR Cluster Survey (MUNICS; 
Drory et al. 2001). The sample includes all the galaxies 
redder than R-K’=5 corresponding to a surface density of 
EROs of 0.04 arcmin“^ at K’< 18. In Tab. 1 the relevant 
photometric information for the two targets is summa¬ 
rized. The two EROs are placed in a high Galactic lati¬ 
tude field selected to contain no bright galaxies or known 
nearby clusters of galaxies. We visually inspected the K- 
band image without finding any candidate leasing objects. 
Thus, we are confident that the magnitudes of the two 
EROs are not affected by magnification of foreground mas¬ 
sive objects. 

Spectroscopic observations, carried out during the 
night 6 March 2002 with 1.5" wide slit, are based on the 
Amici prism dispersing element (Oliva 2001) mounted at 


the near-IR camera NICS of the Italian 3.6 m Telescopio 
Nazionale Galileo (TNG). The Amici prism provides the 
spectrum from 0.85 /rm to 2.4 /rm in one shot with a nearly 
constant spectral resolution of AA/A ~ 50 (l"slit) over 
the whole spectral range. This results in a dispersion of 
~30A (100 a ) per pixel and a full-width at half-maximum 
(FWHM) of ~ 200A (~ 400A) at lOOOOA (20000A). This 
very low-resolution mode is best suited to detect contin¬ 
uum breaks resulting very efficient in identifying old stellar 
systems at z > 1.2. The targets were acquired by means of 
a nearby brighter reference object put in the slit. Dithering 
of the targets along the slit in a A-B-B-A pattern with 
small offset about each of the two positions were used. 
Integrations of 2 minutes for each exposure were adopted 
for all the observations in order to be background limited 
in K’ band. A total integration time of ~100 minutes for 
each source was accumulated. 

After the sky subtraction, the frames have been aligned 
and co-added. Wavelength calibration was performed us¬ 
ing the deep telluric absorption features. The telluric ab¬ 
sorptions were then removed by dividing each of the object 
spectra by an AO reference star spectrum taken at similar 
airmass and adjacent in time. The intrinsic features and 
shape of the reference star were then removed by mul¬ 
tiplying the spectrum by a synthetic AO star spectrum 
smoothed to the appropriate resolution. 

The final extracted spectrum of S7F5-45 and of S7F5- 
254 is shown in Fig. 1 (upper panel and lower panel re¬ 
spectively) along with the relevant mean sky residuals. 
Filled symbols are the photometric data in the I, J and 
K’ bands (see Tab. 1) from the MUNIGS survey (Drory et 
al. 2001). The spectra are normalized to the K-band flux. 
As an independent check on the near-IR photometry, we 
derived a flux calibration of the spectra by comparing the 
mean count rate at K in the targets to the mean count 
rate of the AO star whose K magnitude is known (Hunt et 
al. 2000). For both the EROs the flux density we inferred 
from the spectra matches the near-IR photometry of the 
MUNIGS catalog within a factor of 1.3 in flux. 

3. Analysis of the Amici spectra 

In the following sections we analyze the Amici spectral 
data. We derive the redshift range of the two galaxies 
from their continuum shape and an estimate of their stel¬ 
lar mass based on their K-band luminosity and on the 
mean local mass-to-light ratio. 

3.1. Redshift estimate 

The spectrum of both the EROs is almost flat at wave¬ 
length A > 1.0 /rm and drops rapidly below this wave¬ 
length. We do not detect any emission feature and we es¬ 
timate a limiting equivalent-width (EW) for the emission 
line detection of ~125A and ~150A in the J band for 
S7F5_45 and S7F5_254 respectively. The only continuum 
feature is the change of the slope at A < 1.0 iiui where 
the continuum flux drops by a factor of 4 at A ~ 0.9 /rm. 
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Table 1. Photometry of the two EROs S7F5_45 and S7F5_254 from the MUNICS catalog. All magnitudes are in the 
Vega system and measured within 5” diameter aperture (Drory et al. 2001). 


Object 

RA(J2000) 

Dec(J2000) 

V 

R 

1 

J 

K 

S7F5.254 

13 34 59.6 

16 49 10.7 

25.4±1.2 

24.4±0.5 

23.1±0.7 

19.8±0.1 

17.8±0.2 

S7F5.45 

13 34 25.0 

16 45 48.6 

24.2T0.4 

23.5±0.3 

22.2T0.3 

19.6±0.1 

17.65±0.08 



1 1.5 2 


Fig. 1. Upper panel: NICS-Amici low resolution spectrum 
of the ERO S7f5_45 (thin histogram). The thick histogram 
is the spectrum heavily smoothed to show the continuum. 
The spectrum is normalized to the flux in the K-band. The 
shaded areas represent the atmospheric windows charac¬ 
terized by an opacity larger than 80%. The filled symbols 
are the photometric data in the I, J and K’ bands from 
the MUNICS survey (Drory et al. 2001). The thin line is 
the optical spectrum of Arp220 redshifted at z=1.2 and 
normalized to the J-band flux. The lower histogram is the 
mean sky residuals (with an offset -1.5T0“^® erg cm“^ s“^ 
A“^) extracted above and below of the target spectrum. 
Lower panel: NICS-Amici low resolution spectrum of the 
ERO S7f5_254. Symbols are as in the upper panel. 


We rule out the possibility that the break is due to the 
Lyman-limit because of the detection of both the galaxies 
in the V and R bands. It is worth noting that the broad 
band colors of S7F5_45 could be still compatible, even if 
marginally, with a high redshift (2 ~ 6) galaxy (for in¬ 
stance an SB2 starburst galaxy from Kinney et al. 1996), 
given the large errors in the optical magnitudes. However, 
in this case, the shape of the near-IR spectrum rules out 
this hypothesis being it not compatible with that of any 
galaxy at such high redshift. Dust reddening produces 
smooth spectra (e.g. Schmitt et al. 1997) as it is almost in¬ 
versely proportional to the wavelength: there is no way to 


produce a sharp break in the continuum like the one seen 
at 0.9 < A < 1.0 /rm in our spectra (see also §4). Thus, we 
identify the detected break as the 4000A break produced 
by the G and K stars. This places the two EROs at red¬ 
shift 1.2 < 2 < 1.5 and defines them as possibly early-type 
galaxies. It is worth noting that using the color-color diag¬ 
nostic diagram suggested by Pozzetti & Mannucci (2000) 
both the EROs would lie on the elliptical side of the (J- 
K)-(R-K) plane. For comparison, in Fig. 1 (thin line) the 
optical spectrum of Arp220 (UZC; Falco et al. 1999) red- 
shifted at z=1.2 is superimposed to the Amici spectra. 

The spectra shown in Fig. 1 are comparatively differ¬ 
ent. The drop in the observed flux is sharper in S7F5_45 
than in S7F5_254 and it is well constrained at 0.9 < A < 
1.0 /im. On the contrary S7F5_254 shows a first smooth 
decrease of the flux at<1.0<A<1.3 ^m followed by the 
drop which extends at A < 0.9 /rm. These features suggest 
a redshift of S7F5_254 slightly lower than that of S7F5_45 
and are typical of old stellar populations. 

3.2. Luminosity and stellar mass estimate 

In order to obtain a lower limit to the luminosity of the 
two galaxies and a model independent estimate of their 
stellar mass, we conservatively placed both of them at 
2 = 1.2. Using a k-correction AK=—0.6 mag (Mannucci 
et al. 2001) we estimate a lower limit to their rest-frame 
K-band absolute magnitude of M^ ~ —26.6 mag. Cole et 
al. (2001) find M’^ = —24.4 mag (scaled to Ho=65 Km 
s“^ Mpc“^) for the local luminosity function of galaxies. 
Thus, the luminosity of the two EROs is L ~ 7L*. The 
brightest galaxies in the sample of Cimatti et al. (1999) 
have luminosities lower than L* while the candidate old 
galaxy Cl 0939-I-4713B revealed by Soifer et al. (1999) has 
a luminosity comparable to those of our targets if the pos¬ 
sible magnification induced by the foreground cluster is 
neglected. Thus, the two EROs presented here are among 
the most luminous and, possibly, the most massive evolved 
galaxies detected so far of 2 > 1. If we assume the local 
mean stellar mass-to-light ratio obtained by Cole et al. 
(2001) with the Kennicut IMF (O. 73 M 0 /L 0 ) we derive a 
stellar mass for the two EROs Aistars — 7x 10^^ M 0 which 
would exceed 10^^ M 0 if we assume the Mstars/L they 
derived with the Salpeter IMF (I. 32 M 0 /L 0 ). However, 
these mean mass-to-light ratios are relevant to local galax¬ 
ies. Higher redshift galaxies, on average younger than the 
local ones, should be described by lower ratios. This is 
quantitatively shown in Fig. 2 where the K-band stellar 
mass-to-light ratio as a function of age is shown for dif¬ 
ferent initial mass functions (IMF; upper panel) and for 
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Table 2. Parameters used to construct the grid of models 


SFH T [Gyr] 

0.1, 1, 2, 3, 5, 15, SSP, cst 

Metallicity 

0.2 Zq, 0.4 Zq, Zq, 2 Zq 

A_V [mag] 

04-6 (step of 0.25 mag) 

Extinction law 

Calzetti et al. (2000) 


different star formation histories and metallicity (lower 
panel). 

4. Comparison with population synthesis models 

In this section we analyze the two galaxies in more detail 
adding to the spectral information the available broad¬ 
band optical and near-IR photometry. The aim of this 
analysis is to obtain a more accurate estimate of the red- 
shift of the two EROs and to set a lower limit to their 
stellar mass. To this end we constructed a grid of tem¬ 
plates based on the population synthesis models of Chariot 
& Longhetti (2001). The grid constructed takes into ac¬ 
count different star formation histories (SFH), metallicity 
and extinction. The parameters of the models are sum¬ 
marized in Table 2. The star formation histories consid¬ 
ered, beside the simple stellar population (SSP) and the 
constant star formation (cst), are described by an expo¬ 
nentially declining SFR with an e-folding time r. For each 
star formation history a set of 51 synthetic templates in 
the range 10®-1.5T0^° yr was generated with a Salpeter 
IMF (O.IM© < Ai < IOOM 0 ). By varying the metallicity, 
the star formation history and the extinction we built-up 
a grid of about 38000 templates. 

The redshift of the two EROs was formally measured 
by fitting the drop in the observed low resolution near- 
IR spectra beside the available optical and near-IR pho¬ 
tometric data. The drop in the observed spectra has 
been described by four additional photometric points de¬ 
rived from the observed spectra in the wavelength range 
0.9 < A < 1.2 /im. Another additional photometric point 
has been derived in the H-band. Practically, we convolved 
the flux calibrated spectra with four rectangular 0.1 /rm 
wide filters centered at Ai = 0.95 /rm, A 2 = 1.0 /rm, 
A 3 = 1.08 /xm, A 4 = 1.15 /xm respectively and one 0.2 
/xm wide filter centered at A 5 = 1.65 /xm. The absolute 
calibration of the filters were derived by convolving them 
with the spectrum of an AO star. We used the software 
hyperz (Bolzonella et al. 2000) to obtain the best fit to 
the 10 photometric points for each of these models. 

S7F5_254 - The best-fit value to the redshift of S7F5_254 
is Zbest = l- 22 i(Q' 2 ^ with a minimum reduced = 0.35 
(P(x^)~ 0.96). The quoted errors are not the uncertainties 
in the fit. They represent the boundary of the redshift 
range spanned by the statistically acceptable fit (P(x^)> 
0.68) obtained with the different models. This is shown in 
Fig. 3 where the best-fit value to the redshift as a function 
of the relevant is plotted for those models providing a 
statistically acceptable fit. 


Formally, the best-fit to S7F5_254 data is given both 
by a 10 Gyr old SSP and a 10 Gyr old r = 0.1 model 
with Z=Zq and Ay=0 (see Tab. 3). All the solar and sub¬ 
solar metallicity models with t < 3 provide a good fit 
(P(x^)> 0.9) invoking very old stellar population (> 10 
Gyr) while none of them provide an acceptable fit with 
younger ages. Of course, this does not make sense since 
these ages are much larger than the Hubble time at 
z ~ 1.2. By forcing galaxies to have ages lower than the 
Hubble time at this z the only good result is obtained 
with Z=2Zq models which provide a good fit to the data 
in the case of a 3.5 Gyr old r = 0.1 (P(x^)— 0.95) with an 
extinction Ay=0. This is the youngest mean age we were 
able to obtain. It is worth noting that these properties are 
also displayed by the reddest early type cluster members 
at z ~ 1.2 (Rosati et al. 2000). The best fitting model 
parameters are summarized in Table 3 together with the 
2nd best set of results. In Fig. 4 the best-fitting model, 
the observed spectrum and the photometric data relevant 
to S7F5_254 are shown (upper panel). For comparison, in 
the lower panel, we have superimposed to the data the 
mean observed spectrum of local ellipticals by Mannucci 
et al. (2001, ManOl). The remarkable agreement is indica¬ 
tive once again for the old age of this galaxy. To derive a 
lower limit to the stellar mass we consider the youngest 
population model consistent with the data, i.e. the 3.5 Gyr 
T = 0.1 model. The lowest M/L ratio describing a stellar 
population 3.5 Gyr old fitting the data is ~ 0.4 Mq/Lq, 
obtained with the Scalo and Miller-Scalo IMFs. Since the 
K-band absolute magnitude of S7F5_254 resulting from 
the fitting models is M^k = —27.08 we derive a minimum 
stellar mass = 6 x 10^^ Mq, where we have used 

= 3.4 (Allen 1973). Being the r = 0.1 the youngest 
model, it requires the lowest formation redshift, z/, to fit 
the data. Indeed, any other exponentially decaying burst 
model would require star formation to have begun at a 
larger redshift. Thus, given the redshift of S7F5_254, the 
youngest best-fitting model requires for this galaxy a for¬ 
mation redshift z/>3.5. 

S7F5_45 - The best-fit value to the redshift of S7F5_45 is 
Zbest = 1.46 ± 0.02 (x^ = 0.33, P(x^)~ 0.97). Also in this 
case this value is matched by all the models providing a 
statistically acceptable fit (P(x^)> 0 . 68 ) as shown in Fig. 
3. In agreement with the results obtained in §3.1, S7F5_45 
has a redshift higher than S7F5_254. 

For this galaxy all the exponentially t > 0.1 de¬ 
caying models do not provide acceptable fit to the data 
(P(X^)<< 0.68). The single-burst models provide good fit 
independently on the assumed metallicity. The best-fitting 
model is the Z=O. 4 Z 0 single-burst, 0.5 Gyr old and with 
an extinction Ay=1.5 mag (see Table 3). All the other 
SSP with different metallicities provide good fit as well 
(0.93<P(x^)< 0.96) giving mean ages in the range 0.25 
Gyr {Z=2Zq) and 0.7 Gyr (Z=O. 2 Z 0 ) and an extinction 
1.25<Ay < 1.5. Only the t = 0.1 {Z=Zq) is the exponen¬ 
tially decaying model providing a good fit (P(x^)— 0.95). 
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Table 3. Parameters resulting from the best fits of stellar population models to the near-IR spectra and photometric 
data 


ERO ID 

SFH 


P(x") 

Z 

Age 

(Gyr) 

Av 

(mag) 

Z 

(Z©) 

Mk 

(mag) 

/L 

(Mo/Lo) 

'■^^stars 

(IO^Mq) 

Best Fit 
S7F5_254 

r = 0.1,SSP 

0.35 

0.96 

1.22 

10 

0 

1 

-27.05 

1.26 

19 

S7F5_45 

SSP 

0.33 

0.97 

1.46 

0.5 

1.5 

0.4 

-27.69 

0.4 

10 

2nd Choice 
S7F5_254 

T = 0.1 

0.36 

0.95 

1.22 

3.5 

0 

2 

-27.08 

0.46 

7 

S7F5_45 

r = 0.1 

0.37 

0.95 

1.45 

0.5 

1.75 

1 

-27.73 

0.3 

8 


z 

30 20 10 5 2 1 0.5 0 

1 .5 

1,4 

1 .3 

1.2 


1,1 


1,4 

1 .3 

1.2 


1 .1 



108 109 IQIO 

age [yr] 



0.2 0,4 0,6 0,6 


(P(y2)S0.68) 


Fig. 2. Upper panel: Mass-to-light ratio as a function of 
age for different IMF and a star formation history de¬ 
scribed by a SSP. Lower panel: Mass-to-light ratio as a 
function of age for Salpeter IMF, different metallicity and 
a star formation history described by SSP (upper curves) 
and by cst (lower curves). 


In this case we obtain a mean age of 0.5 Gyr and an ex¬ 
tinction Ay=1.75 mag. In Fig. 5 the best-fitting model, 
the observed spectrum and the photometric data relevant 
to S7F5_45 are shown (upper panel). These results sug¬ 
gest that the emission of this galaxy is dominated by a 
population of stars recently formed. On the other hand 
this population has to be formed in a burst with a time- 
scale much shorter than the mean age itself (i.e. t ^ 0.5 
Gyr) otherwise we should obtain an acceptable fit also 
with T > 1 models. Unless to hypothesize that all the 
stellar mass of this galaxy is formed in such a short burst, 
the population of stars recently formed has to be superim¬ 
posed to an older population. This is quantitatively shown 
in Fig. 5 (lower panel) where a composite stellar popula¬ 
tion (GSP) is superimposed to the data. The GSP shown 
is the weighted sum of the best-fitting SSP and an old 
(5 Gyr, r = 1, Z=Zq, Ay=0) template spectrum. They 


Fig. 3. Best-fit value to the redshift of S7F5_254 (upper 
panel) and of S7F5_45 (lower panel) as a function of the 
relevant reduced for slf th® models providing a statis¬ 
tically acceptable fit (P(x^)> 0.68). 

have been weighted so that the old spectrum contribute 
for 60% of the total stellar mass. It is evident that the 
population of stars recently formed dominates the emis¬ 
sion even in this extreme case. Thus, the stellar mass-to- 
light ratio relevant to a mean age of 0.5 Gyr represents 
the most “conservative” choice we can make. The low¬ 
est ratio for this mean age is ~ 0.15 Mq/Lq, the value 
given by the Miller-Scalo IMFs (e.g. Fig. 2). Given the 
K-band absolute magnitude of S7F5_45 resulting from the 
fitting models (M^y = —27.7) we derive a minimum stellar 
mass = 4 x 10^^ M©. For this galaxy, the best¬ 

fitting model requires a lower limit to the formation red- 
shift Zf>2. 

5. Summary and Conclusions 

We have presented low resolution near-IR spectra of two 
candidates z > 1 massive {Mstars > 10 ^^ Mq) galaxies 
out of the 15 candidates selected from the MUNIGS survey 
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^obs [Mm] 


Fig. 4. Upper panel: The best fitting template (thick grey 
line) is superimposed on the observed smoothed spectrum 
(thin black histogram) of S7F5_254. The filled circles are 
the photometric data in the V, R, I, J and K’ bands from 
the MUNICS survey. The filled triangles are the additional 
photometric points derived by the observed spectrum (see 
§4 for details). Lower panel: The mean observed spectrum 
of local ellipticals (Mannucci et al. 2001, ManOl) is super¬ 
imposed on the observed smoothed spectrum of S7F5_254. 
Symbols are as in the upper panel. 


(R-K>5 and K< 18). The spectra, obtained with NICS- 
Amici, place the two EROs at 1.2< 2 < 1.5. The stellar 
mass derived by their rest-frame 2 = 1.2 K-band absolute 
magnitude (M/f ~ —26.6) and assuming the local stel¬ 
lar mass-to-light ratio (O.73 MqL0 ) is Mstars ^ 7 x 10^^ 
Mq. In order to better constrain the redshift and to ob¬ 
tain a lower limit to the stellar mass we have compared 
the spectra and the photometric data with a set of tem¬ 
plate spectra by means of a best-fitting procedure. 
The best-fit to the data provides a redshift 2 ~ 1.22 
for S7F5_254 and 2 1.46 for S7F5_45. We estimate a 

lower limit to their stellar mass of = 6 x 10 ^^ Mq 

and MTtTrs = 4 X 10^^ Mq respectively, under the most 
conservative results provided by the best-fitting models. 
Thus, these two galaxies are among the most luminous 
and massive evolved galaxies detected so far at redshift 
2 > 1. The galaxy S7F5_254 is comparatively older than 
S7F5_45. The youngest mean ages of the population of 
stars dominating the emission of the two EROs are 3.5 Gyr 
and 0.5 Gyr respectively which suggest formation redshift 
of 2 /> 3.5 and Zf>2. These results, even if limited to two 
EROs, point to a galaxy formation scenario where massive 
evolved galaxies are already well in place at 2 > 1 . 5 . Given 
the high values of the lower limit to the Mstars of our 
EROs, they could severely constrain the galaxy formation 
models. Spectrophotometric models of galaxies with an 
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Fig. 5. Upper panel: The best fitting template (thick grey 
line) is superimposed on the observed smoothed spectrum 
(thin black histogram) of S7F5_45. Symbols are as in Fig. 
3. Lower panel: A composite stellar population (CSP) is 
superimposed to the data. The GSP is the weighted sum 
of the best-fitting SSP and an old (5 Gyr, r = 1, Z=Z 0 , 
Av=0) template spectra. They have been weighted so 
that the old component contribute 60% of the total stellar 
mass. 

accurate modeling of the horizontal branch and includ¬ 
ing different values of a/Fe (Maraston & Thomas 2000; 
Thomas, Maraston & Bender 2002) as well as the emis¬ 
sion from the gas component (Gharlot & Longhetti 2001), 
will allow us to interpret the higher resolution near-IR 
spectroscopic data we expect in the fall next year at VLT- 
ISAAG. These data make it possible to restrict the ranges 
of ages and metallicity providing a more accurate estimate 
of the stellar mass of the two massive galaxies presented 
here. Moreover, our XMM observations scheduled in 2003 
will allow us to put severe constraints on the presence of 
AGNs in early-type galaxies at high redshifts. 
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